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This paper is a review of the present state of knowledge of magneto-telluric fields. The 
subject has to do with the combined analyses of the geomagnetic and the telluric (earth- 
current) fields on the surface of the earth. Usually, the objective of such investigations is to 
obtain information about the earth's crustal layers. However, for a sensible use of the 
method it is desirable to understand something about the source of the fields. In this paper, 
the various suggestions for the source mechanism are discussed. Then a fairly detailed 
review of previous work on the theory of the magneto-telluric interpretation is given. In- 
cluded are a number of three-layer interpretation curves. The influence of earth curvature 
in magneto-telluric theory is treated in a mathematical appendix which is, in itself, a self- 
contained derivation of the various formulas. 

1. Introduction 

The temporal variations of the geomagnetic field have been investigated extensively for 
many years. However, only quite recently have the rapid magnetic variations been studied in 
conjunction with the variations of the telluric (earth-current) fields. In fact, only some twelve 
years ago, Tikhonov [1950] in the USSR, and Kato and Kikuchi [1950] in Japan, pointed out 
that the electrical characteristics of the deep strata of the earth's crust could be determined 
from a combined analysis of geomagnetic and telluric field variations. Since then, a large 
number of related investigations have been carried out, particularly in the USSR. It is the 
purpose of this paper to present a self-contained account of the theory of the phenomena. 
Sufficient curves are given to provide a basis for interpretation of experimental data. 

2. Nature of the Sources 

The actual mechanism which produces the short-period variations (i.e., frequencies of the 
order of 1 to 10~ 3 c/s) is not yet well understood. In recent years it has been suggested on 
numerous occasions that the phenomena are related to magneto-hydrodynamic (i.e., MHD) 
waves in the exosphere or the ionosphere of the earth. 

Holmberg [1951], Lehnert [1956], and Maple [1959] have proposed that the sources of the 
observed micropulsations arise from intralayer MHD oscillations within the E and F layers of 
the ionosphere. The principal objection to such a hypothesis is the high attenuation of MHD 
waves in such regions as pointed out by Piddington [1959] and Ellis [I960]. 

Dungey [1955] has suggested that MHD waves form a standing wave pattern along the 
geomagnetic field lines. In the fundamental mode the period of oscillations would then be 
expected to increase with geomagnetic latitude. Some support for this view has been given 
by Obayashi and Jacobs [1958], who found that the period of some micropulsations did appear 
to increase with latitude. However, Ellis [1960], in Australia, found that there was no 
observable change in the period for three stations ranging from 28° S to 51° S geomagnetic 
latitude. Dungey [1955] has also suggested that instabilities of the Chapman-Ferraro outer 
boundary of the geomagnetic field could generate waves which propagate downward as MHD 
waves. It would be expected that such waves would tend to reach the ionosphere in auroral 
regions and consequently micropulsation from this source would not be observed at middle 
and equatorial latitudes. However, they might be propagated in horizontal MHD ducts as 
envisaged by Bomke et al. [I960]. 
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Quite recently, Gallet [1959] and Ginzburg [1961] have suggested that solar corpuscular 
streams produce magneto-hydrodynamic waves. In particular, Ginzburg considers the radia- 
tion emitted by an isolated ion moving with velocity u along a magnetic field H in a magneto- 
active plasma. Two cases are distinguished, (1) slower-than-light motion, when the ion 
velocity u is less than the phase velocity v Ph of the emitted wave, and (2) f aster- than-light 
motion, when the ion velocity is greater than the phase velocity of the emitted wave. 

In the first case, mentioned above, the radiation results from the Larmor precession of 
the ion around the lines of force of the earth's steady magnetic field H. The angular fre- 
quency of the rotation is Q,i=eHjMc, where e is the ion charge, M is the ion mass, and c is 
the velocity of light in vacuum. For a fixed observer in the plasma, the frequency of this 
radiation is 

w 1- (u/c)N cos d W 

where N is the refractive index of the plasma at frequency a/, and 6 is the angle between the 
direction of wave propagation and the ion velocity u. This is really the normal Doppler 
effect and it results from the motion of the emitter which, in the present situation, is the 
gyrating ion. 

In the second case, denoted (2) above, the motion is faster-than-light and the frequency 
of the emitted radiation is given by 

0< 



-Ncosd- 
c 



(2) 



This has been called the anomalous Doppler effect [Frank, 1942]. Contrary to the normal 
Doppler effect, in the latter case a/ may be less than Q t if u is sufficiently great. 

By using known expressions for the refractive index of a plasma, and specifying the 
direction of wave propagation, the frequencies a/ emitted by the ion can be calculated. To 
simplify the discussion the waves are assumed to propagate along the earth's magnetic field 
(i.e., longitudinal propagation). Then, on the further assumption that the collision frequencies 
are negligible [Fejer, 1960; Wait, 1961], 

iV 2 ^l ^ ^ - (S) 

where oj ,i and co , e are the ion and electron plasma frequencies, respectively. 
In the earth's exosphere and for co<<^ 6 , the magnitude of Nis large. Then 

N^ 7 ^ ,- where K =^- (4) 

, , ratio of electron mass 1 » , , . mu £ 

and where a=— —. ^ =77^77 for hydrogen ions. I he last magneto-sonic wave 

ratio of ion mass 1836 J & 

corresponds to the adoption of the upper sign in the above expression. The Alfven wave 

is obtained when the lower signs are used. It may be noted that for zero frequency the 

velocity v A of the Alfven wave is given by 



nV < 6 > 



Inserting the equation for N into eqs (1) or (2), one easily finds that 

<'*2)[«>(>-£)7 
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where the upper sign is to be chosen for slower-than-light motion and the lower sign for faster- 
than-light motion. 

Ginzburg [1961] has shown that the intensity of the cyclotron radiation is identically 
zero for the slower-than-light case when the propagation is longitudinal (i.e., 6=0). Thus, 
the only case of physical interest is when the lower or + sign is used in the above equation. 

Actually, a charge moving along a helical path also radiates harmonics of the Larmor 
frequency 12*. Thus, on replacing 12* in eq (6) by pQ t where p is an integer, we find that the 
equation, for faster-than-light motion, becomes 

This equation may be used to obtain the frequencies of emission from the specified value 
of u/v A . 

A specific example is now considered. At a distance of several earth radii H is of the 
order of 5X 10~ 3 gauss * and if u^7X 10 6 in/sec, n^ 107cm 3 , then the Alfven velocity v A ^ 10 6 
m/sec and 12^50 radians/sec. Solution of eq (7), when p=l, leads to the values co/2x=0.77, 
777, and 13,000 c/s. It is only the lowest frequency that would be of interest in magneto- 
telluric investigations. 

In general, when the conditions for longitudinal propagation no longer hold, the situation 
is much more complicated. In this case Ginzburg [1961] has shown that the radiation sped nun 
extends over a much wider range. He estimates that the excitation zone extends from about 
13 earth radii down to about 1.5 earth radii. The lower limit depends on collision damping, 
noncollision cyclotron absorption, and the increase of the Alfven velocity to unfavorably 
large values. 

Magneto-hydrodynamic waves of the type described above are probably the main source 
of the short-period geomagnetic field variations. On the basis of Ginzburg's analysis one 
would expect the spectrum to peak in the region of 1 c/s or so. However, observations by 
Campbell [1959], and Berdichevsky and Brunelli [1959] on the surface of the ground indicate 
that the spectrum peaks around 0.03 c/s with a very rapid decrease for frequencies above 
about 0.2 c/s. The apparent discrepancy can be reconciled if the transmission characteristics 
of the ionosphere are considered [Akasofu, I960]. 

For longitudinal propagation within the ionosphere the complex propagation constant T is 
given by Hines [1953] and Dungey [1955] 

r.2~ ^ -ivp/a 



(HWp,) (lTiv e /tt e )(l±ii> p /tt p ) 



where p p is the mass density of the ions, v e is the effective collision frequency from the electrons 
with the neutral particles, and v v is the effective collision frequency for the positive ions. 
The transmission coefficient T for the ionosphere is then given approximately by 






(9) 



where h x and h 2 are the effective heights of the bottom and top of the ionosphere. This formula 
gives the ratio of the field below for a wave of unit amplitude incident from above. It is only 
valid when the medium is slowly varying. For frequencies 0.01, 0.1, and 1 c/s, the corre- 
sponding ratios are 0.32, 2.8X10 -2 , and 8.2X10 -6 if Dungey's model of the ionosphere is used. 
This indicates that the ionosphere is almost transparent for frequencies less than about 0.1 
c/s. For the higher frequencies, the ionosphere becomes almost opaque. In fact, if the main 
source of micropulsations is above the F layer it means that the sources must be very intense. 



*The term gauss is actually a magnetic flux density, although it is frequently used in the literature for magnetic field. Note that for magnetic 
flux, one weber/m2— 10* gauss -10 3 gammas. 
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For these "higher" frequencies it is more probable that the sources are currents at E region 
heights in the ionosphere produced by some other mechanism [Baker and Martyn, 1953; 
Baker, 1953]. 1 Also, at frequencies of the order of 1 c/s and higher, one may expect contri- 
butions from lightning discharges which are known as "sf erics." Support for this view is 
obtained from the work of Holzer and Deal [1956], who find that there is a strong correlation 
between worldwide thunderstorm activity and the electromagnetic signals in audio and sub- 
audio frequency range. 

Because of the uncertainty and the conjectural theoretical nature of the source mecha- 
nisms, one could feel very insecure in attempting to interpret field measurements on the earth's 
surface in terms of crustal features. However, if the ratios of the mutually orthogonal com- 
ponents of the tangential electric and magnetic fields are considered, some very definite con- 
clusions can be drawn about the crustal structure of the earth. Such studies are usually referred 
to as magneto-tellurics, since the analyses of geomagnetic and telluric (earth-current) variations 
are combined. 

3. Review and Discussion of Previous Work on Magneto-Tellurics 

Most investigations of magneto-telluric fields boil down to a study of the interrelation 
between the tangential components of the horizontal electric and magnetic fields at the surface 
of the earth. As far as this writer is able to ascertain the first definitive paper dealing with 
this subject appeared in 1950, and was authored by Tikhonov [1950]. He realized that rapid 
geomagnetic variations and earth currents, observed at the surface of the earth, must be con- 
nected by some definite relationship. He showed that, at low frequencies, the amplitude of 
the derivative of the component H x of the magnetic variations is proportional to the (orthogo- 
nal) component of the electric field E y . 2 This was in agreement with the experimentally 
established fact that there is a proportionality between these quantities. Tikhonov's model 
of the earth's crust is a planar layer O^z^l of finite conductivity a lying upon an ideally con- 
ducting substrate. Implicitly in his analysis, it was assumed that horizontal gradients of the 
fields could be neglected. Thus, for a spectral component of frequency co he found that, at 

2=0, 

in uH x g*Eyy coth (7/) (10) 

where 

y= (icrnow)* and /z =4X 10" 7 . 

Using previously published data on the observed diurnal variations at Tucson (Arizona) and 
Zui (USSR), Tikhonov computed the value of a and I which best filled the first four har- 
monics. For Tucson, the values were about 4X10" 3 mhos/m and 1,000 km, respectively. For 
Zui, the corresponding values were about 3X10" 1 mhos/m and 100 km. 

In a later paper [Tikhonov and Lipskaya, 1952], the horizontal gradient of the fields was 
allowed to be finite (although the earth's crust was still regarded as a horizontally stratified 
medium). The authors postulated that the field components of long period may be represented 
as a wave which is propagated from east to w^est with the velocity of the earth's rotation. 
Using the same data as mentioned above, they obtained revised estimates for the conductivity 
a of the upper stratum of thickness I under the assumption of an ideally conducting substrate. 
For Tucson, the values were approximately 10~ 2 mhos/m and 1,100 km, where for Zui the cor- 
responding values were about 7X10 -1 mhos/m and 110 km. In this paper they also showed 
that the measured values of the vertical magnetic field variations were consistent with the model 
and the postulation of linear east-to-west motion of the fields. 

In the third paper of this sequence [Lipskaya, 1953], the electromagnetic equations for the 
model described above were cast in a form to clearly demonstrate the various relationships 
between the field components. 



1 W.H.Campbell has recently proposed that small, short-period changes in the gross ionospheric pattern within the auroral electrojet are a 
principal cause of micropulsations (private communication, May 1962). 

2 The Cartesian coordinate system is chosen so that the z axis points downward and the earth's surface is z=0. 
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Subsequent to the Russian work mentioned above, Cagniard [1953] published a paper 
which has been extensively referenced since. His analysis which assumed plane wave incidence, 
developed formulas which related H x and E y on the surface of a stratified conducting medium. 
A discussion of the limitations in Cagniard 's results appeared shortly thereafter [Wait, 1954]. 
The essential point made by this writer is that the proportionality between H x and E y is only 
valid if the fields themselves do not vary appreciably in a horizontal distance of the order of a 
"skin depth" in the ground. In the case of a homogeneous flat earth it was indicated that this 
distance was of the order of |7 -1 |. For example, at a frequency of 10~ 3 c/s, a ~10 -3 mhos/m, 
|7 _1 | is about 350 km. Consequently, the field should be uniform over a considerably broad 
area to permit the Cagniard interpretive procedure to be applied. This limitation is, of course, 
in addition to the requirement that the crust al layers themselves are uniform. Very recently 
Price [1962] has indicated that the limitation mentioned [Wait, 1954] becomes much more 
stringent when the magneto-telluric method is applied to a stratified earth in certain important 
instances. 

In a later paper, Tikhonov and Shakhsuvarov [1956] discussed the methods for calculating 
the admittance H x /E y at the surface of a horizontally stratified earth of any number of layers. 
Again, horizontal field gradients were neglected. Actually, the formulas given were a special 
case of an earlier general analysis [Wait, 1953a and b] where the incidence was oblique. 
Tikhonov and Shakhsuvarov [1956] give a few curves of the ratio H x /E y for both two- and 
three-layer structures. Some asymptotic approximations of the impedance formula for parallel 
layers were discussed by Berdichevsky and Brunelli [1959]. The derivation of the basic formu- 
las is very similar to that found in earlier papers [e.g., Wait, 1953a]. 

The problem has also been investigated by Scholte and Veldkamp [1955] in Holland. Their 
analysis of the variations of H x and E y is essentially the same as those of Tikhonov, Cagniard, 
and others. By a relatively simple method of data analysis they estimate the amplitude and 
phase of H r jE y without actually performing spectral analyses. They used the experimental 
data from the magnetic observatory at Witteveen for frequencies in the range from lO -4 to 
10" 1 sec. The resulting curves of amplitude and phase versus frequency fitted a two-layer 
earth model with the upper conductivity <j y =2 mhos/m and the lower conductivity o^^lO -1 
mhos/m. The thickness of the upper stratum was 6^ = 600 m. In the same paper Scholte and 
Veldkamp also discuss the influence of earth curvature on the ratio H V IE X . They conclude 
that the effect is very small. (A more general theory is given in the appendix.) 

An analysis of the continuously varying conductivity profile was carried out by Bossy 
and De Vuyst [1959] in Belgium. They showed that for normal incidence the ratio II x IE y on 
the surface could be expressed in closed form when the conductivity <r(z) varied with depth in 
the manner 

a(z) = a (l+^~" (11) 

where tr and a are constants. They used such a model to interpret experimental data of the 
phase of E x /H y for frequencies from about 3 X 10~ 3 to 1 c/s at Dourbes. However, they found 
it necessary to modify the model by having a (well conducting) homogeneous surface layer of 
thickness about 500 m overlying the inhomogeneous substrate of poor conductivity (with fi 
about 9). 

In all the work mentioned above, the natural magneto-telluric field has been studied only 
at frequencies less than about 1 c/s. Vladimirov [1960] has described an investigation of the 
use of higher frequencies in the range from 0.3 to 1,000 c/s. Field experiments were made in 
the Ryl'sk district of the Kursk region in the USSR. The area consists of sand and clay 
deposits underlain at a depth of approximately 500 m by a crystalline base of practically 
infinite resistivity. The measured values of the ratios \E y /H x \ 2 and \E X /H y \ 2 on the surface 
were consistent with an assumed two-layer structure where the upper conductivity <ti^5X10~ 2 
mhos/m and the lower conductivity o- 2 ^0. From the experimental curve, the thickness of 
the upper stratum was deduced to be 450 m. In a subsequent paper Vladimirov and Kolmakov 
[1960] discuss the resolving power of the magneto-telluric method. They confine their dis- 
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cussion to a three-layer model with an infinite resistive basement. The principle of equivalence 
stated by them is that for a given conductivity cr x and thickness hi of the uppermost layer, 
the theoretical three-layer curves of \H X /E y \ practically preserve their shape when the con- 
ductivity a 2 and thickness h 2 of the intermediate layer vary over specified limits. (This ques- 
tion is discussed again in the present paper.) Vladimirov and Kolmakov then conclude that 
the same curve may characterize different geoelectric sections and can be interpreted only 
with known values of a 2 and h 2 . Some further comments on the same subject are given by 
Vladimirov and Nikiforova [1961] and Vladimirov and An [1961]. 

In an interesting paper, Chetaev [1960], also from the USSR, illustrates a procedure to 
obtain information about the anisotropy of a homogeneous stratum by measuring the surface 
impedance in orthogonal directions. The model he considers is a half -space with a medium 
whose longitudinal conductivity is a t and the transverse conductivity is o> The angle of 
inclination of the anisotropy is a (i.e., a=0 corresponds to a flat lying anisotropic medium 
with a t the horizontally directed conductivity). Chetaev shows that if the electric component 
E y is transverse to the strike of the structure, 

Z=E y /H x =-(i^la t yil + [Mal)-l} sin 2 «]* (12) 

and if the electric component E x is along the strike, 

Z=E x \H y =(i^\<jtf (13) 

which is not affected by <r t Since a t \^i is usually greater than one, the minimum value of 
the surface impedance \Z\, as a function of azimuth, is (/xco/o-^ which occurs when the electric 
vector is along the strike. 

It has been pointed out by Pokityanski [1961] that the interpretation of field results 
should take account of the anisotropy of the structure. In fact, in general, the tangential 
fields are related by 

E x = Z xx H y — Z xy H x (14) 

E y =Z yx Hy—Z y yH x (15) 

where the coefficients are the components of a surface impedance tensor. 

In the case of parallel stratified or flat-lying media Z xy and Z yx vanish and Z xx =Z yv . Of 
course, the components of the impedance tensor depend on the choice of the (x,y) axes. For 
an anisotropic homogeneous half-space or for an inhomogeneous structure where the con- 
ductivity does not vary along one of the horizontal directions, it is possible to find special 
directions of x and y (denoted u and v, respectively), such that 

E U =Z X H V (16) 

and 

E V =-Z 2 H U (17) 

where Zi and Z 2 are the principal values of the tensor impedance. Consequently, if the u and v 
components of the tangential fields are measured, Z x and Z 2 may be simply calculated. Un- 
fortunately, those principal directions are not always known beforehand. Pokityanski [1961] 
has described an ingenious graphical scheme to determine the u and v directions from measure- 
ments of the nonorthogonality of the tangential E and H vectors. Strictly speaking, the 
method is only valid for source fields which are linearly polarized. A more straightforward 
procedure would be to return to the surface impedance relations and write them in the form 

Z x =EJHy= [Z xx —Z xy a] (18) 

and 

Zy=E y /H X = — [Zyy—Zyja] (19) 

where a=H x /H y and where Z x and Z y are the impedances as measured in the x and y directions, 
respectively. 
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It is clear that Z x and Z y depend on a and, consequently, they will be a function of the source 
field. However, the elements of the surface impedance tensor can be calculated if at least two 
independent measurements are made. Thus, for a particular frequency, these will provide the 
following sets of values, Z XtX , Z yA , a x and Z Xt2 , Z y<2 , and a 2 where the subscript 1 or 2 is used to 
distinguish between the two measurements. Then, provided ai^a 2 , it easily follows that 

Z xy = Z *< l - Z *< 2 , (20) 

OLi— a 2 

z --(i/« 2 )-(i/ ai )' {2l) 

Z xx =Z Xi ! + (Z x ,Jai) 

=Z z , 2 +(ZJa 2 ); (22) 

Zyy — Zy t i -f" {Zyxjai ) 

= -Z Vt2 +(Z yx /a 2 ). (23) 

In principle these equations could be used to determine the elements of the surface impedance 
tensor for an elliptically polarized field. It is probable that the apparent scatter in the mag- 
neto-telluric data, as measured by Garland and Webster [1960] and Watt et al. [1962], is a 
consequence of the anisotropy or of tilting the structure. 

In a recent paper, Kovtun [1961] has discussed, in a general way, the nature of the mag- 
neto-telluric fields for two-dimensional inhomogeneous structures. It is assumed that electrical 
properties vary only in the (x,z) plane; the surface of the earth is the (x,y) plane. It is stated 
incorrectly that the field can be decomposed into an independent set of TM (transverse mag- 
netic) and TE (transverse electric) waves. While this may be true for the incident or primary field, 
it is overlooked that the boundary conditions couple these waves together [Wait, 1959] in the 
general case. Consequently, his subsequent discussion can only be considered approximate. 

Recent activity in the United States on magneto-tellurics is localized mainly at the Massa- 
chusetts Institute of Technology and the University of Texas, and DECO Electronics, Inc., 
in Boulder, Colo. Unfortunately, little of the work has been published in the open literature. 

In a short note, Cantwell and Madden [1960] at M.I.T. report some preliminary magneto- 
telluric measurements in Massachusetts covering the frequency range from 5X10 -3 to 1 c/s. 
The analysis of the records made use of auto- and cross-correlation techniques so that estimates 
of the coherency between electric and magnetic signals were obtained. For the interpretation 
they reject records which have low coherency as they are claimed to result mainly from noise 
due to the instability of the coil mounts. According to Cantwell and Madden [I960], 3 the sig- 
nals associated with high coherency between electric and magnetic signals yield consistent 
estimates of resistivity. They find that their data is compatible with a two-layer model with the 
upper, conductivity o-!^1.2X10~ 4 mhos/m and the lower conductivity o- 2 <1.2X10 _1 mhos/m. 
The thickness h x of the upper stratum is 70 km. 

4. New Numerical Results for Horizontally Stratified Structure 

Although this paper has been primarily a review of the theories and concepts in magneto- 
tellurics, it is considered to be worthwhile to include some specific curves which can be used to 
interpret future experimental data. The numerical results were actually obtained in connec- 
tion with theoretical studies of radio wave propagation in the presence of stratified media. 
Here, for convenience, the results are given in dimensionless form to facilitate and broaden the 
application to field problems. 



3 Apparently the curves in figure 3 of Cantwell's and Madden's paper are mislabeled. The values indicated for p2 should all be divided by 10. 
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Figuke 1. Stratified model of the earth's crust. 

The model is quite simple, the earth is assumed to be horizontally stratified and consists of 
three homogeneous layers. The upper layer is of thickness hi with conductivity a x , the middle 
or intervening layer is of thickness h 2 with conductivity o- 2 , and the bottom layer is of infinite 
thickness and has a conductivity o- 3 . In terms of a Cartesian coordinate system (x,y,z) the 
earth's surface is 2=0 and the interfaces between layers 4 are at z=hi and z=hi+h 2 . The 
magnetic permeability is assumed to be constant throughout and equal to /j . Furthermore, 
displacement currents are neglected in all three conducting layers. The problem is now to 
find an expression for the ratios of the tangential fields E and H at the boundary between free 
space and the earth. If the horizontal gradients of the fields are negligible the result can be 
written (for a time factor e ia}t ) 



ry —Ey E X /^M0<Ai/~>// 7 0-2 0"3 h{\ 



where Q is a function of the four variables indicated. Explicitly, 

A 

_<2+(<r2/Vi)*tanh [(i<rin <a)*hi 
where 



(0-2/0-1)* + $ tanh [(i<7ijUow)*A 2 ] 



Q= 



l + C^M)* tanh [('io- 2 Moa))^2] 
(o- 3 /<r 2 )*+tanh [(i^/iow) 1 ^] 



(24) 



(25) 



(26) 



This result follows directly from a previously derived general formula applicable to any number 
of layers [Wait, 1953a and b]. 



< See figure 1 . 
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Figure 2. Amplitude and phase of the complex quantity Q for a three-layer model plotted as a function of the 

dimensional frequency /depth factor (tri/zco) Ihi. 
[For this set of curves the exciting field is assumed to be uniform (see text) .] 



640G18— 62 2 



517 



o I 




(oj^w) l/2 h, 



TTTT 




J I I II I III I I I I I I III | I I I I I I 



I I Mil 1 I I I I M il I I I I I III 



10" 



rr3 



r2 



I0" 1 



10 



(CTj/iCcP) l/2 h | 

Figure 3. Amplitude and phase of the complex quantity Q for a three-layer model plotted as a function of the 

dimensional frequency /depth factor (<muD)ihi. 
[For this set of curves the exciting field is assumed to be uniform (see text).] 
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Figure 4. Amplitude and phase of the complex quantity Q for a three-layer model plotted as a function of the 

dimensional frequency /depth factor (ai/jico)ihi. 
[For this set of curves the exciting field is assumed to be uniform (see text).] 
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It can be seen from eq (25) that if h { -^™ , Q tends to unity. This suggests that, for the 
general case, one defines an apparent conductivity a a such that 



Thus 



Z=(^)*=(^)V (27) 

Q=(<ril<r a ) i or <rj*i= Q' 2 . 



Clearly, a a must be complex in order to admit such an equivalence with a homogeneous half- 
space. The concept of an equivalent or complex conductivity has been discussed previously 
in connection with radio propagation over stratified media [Wait, 1953b, 1958]. 

Extensive numerical results for Q over a wide range of the four parameters are available 
[Jackson et al., 1962]. These greatly extend the set of curves of \o- a /a r \, for a s =0 only, pub- 
lished by Yungul [1961]. A sample of some of these numerical data is given in this paper in 
graphical form. In the first set shown in figures 2a to 6b, the ratio <t 2 Ig x is fixed at 1/100. 
This corresponds to an intervening layer of very poor conductivity. In each case, the abscissa 
is the dimensionless ratio (o^ocd)^/^. For these five sets of curves, the values of h ] /h 2 are 1/50, 
1/20, 1/10, 1/5, and 1/2, respectively. Thus, in all cases, the upper layer is thin compared with 
the middle or intervening layer. 

It is interesting to note that for the low frequencies [corresponding to small values of 
((Tijitoco)^!], (r a always approaches <r 3 the conductivity of the bottom layer. At the high fre- 
quencies, a a approaches o x as noted above. In the intermediate region an interesting transition 
takes place. The right portion of the curves is roughly characteristic of a two-layer structure 
and the shape is determined mainly by the characteristics of the upper layer. 

In the second set of curves shown in figures 7a to 9b, the ratio h 1 /h 2 is fixed at 1/20 and the 
conductivity ratio a 2 /<ri takes the values 1, 0.1, and 10~ 3 . In the first set shown in figures 7a 
and 7b, the upper two layers combined in a single layer of thickness Qi\-\-h 2 ). For this case, 
the results are identical to the standard two-layer curves [Cagniard, 1953; Wait, 1953a]. The 
interesting thing about the three sets of curves, in figures 8a to 9b, is their similarity. Here 
it appears that the conductivity of the intervening layer a 2 plays a very small role. Actually, 
this is a general characteristic of all three-layer magneto-telluric curves provided a 2 is somew I ml 
less than a { . 

The relative insensitivity of the factor Q to the conductivity and thickness of a poorly 
conducting intervening layer can be demonstrated directly from the basic equations given 
above. For example, under the conditions that 



and thus eq (25) becomes 



(o- 2 /x a;)^A2<<l 

A 

Q^ (0-2/0-3)" 

n^ J + QrsM)* tanh [(wijuo)*Ai] . 

^ = (0-3/a-^+tanh [(WaOtoil ' ^ 



Therefore, under the restriction stated above, the three-layer structure is equivalent to a 
two-layer model whose parameters are a u <r 3 , and hi. The result is essentially independent of 
g 2 and h 2 . 

It is worth pointing out that a dual three-layer model exists which broadens the applica- 
tion of these numerical results. The parameters of the dual model are indicated by primed 
symbols and they are related to the original problem as follows: 

O^/^l — 0"l/02 

0-3/ 'a[= ail V 3 

(<r[iJLou>)Wi= (o-iMo^)^i 

(<T2fXQO>)W 2 = (o-:>Mo^)^2. 
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It then immediately follows from eqs (25) and (26) that 

Q' ( -y/vifJLowhl} -7J -7* T7 )= , — j-r- (29) 

V o-i o"! As/ n / r-7 t o- 2 cr 3 A,A 

\ 0"! <Ti /^2/ 

The apparent conductivity a' a in the dual or transformed problem is given by 

Consequently, the ordinate in figures 2a, 3a, to 9a can be regarded as the quantity \(<r'j a^l 
in the transformed problem. The conductivity ratios are all inverted and 

h[ ci hi 

fl 2 <?2 ^2 

As is obvious, the ordinate in the phase curves in figures 2b, 3b, to 9b is the negative of the 
phase of Q' . 

In the five sets of curves from figures 2a to 6b, the values of h[/h' 2 are thus 2, 5, 10, 20, and 
50, respectively, and the conductivity ratio a^o-J is fixed at 100. 

The dual property of magneto-telluric curves for horizontally stratified structures has 
been discussed recently by Kolmakov [1961]. His analysis seems to be unnecessarily involved. 
Actually, using our correction factors (i.e., amplitude and phase of Q), the duality is almost 
self-evident. 

The curves shown in the preceding figures are based on the important and often over- 
looked assumption that the exciting field is effectively uniform. Stated in another way, the 
fields should not vary appreciably in the horizontal direction. To indicate the significance 
of this assumption for a horizontally stratified medium, it is desirable to consider a two- 
dimensional spatially periodic field. The coordinate system may then be chosen so that b/dy=0 
and d/dx=—ikS where &=27r/free-space wavelength =w/c and S is a dimensionless quantity 
which may be complex. Therefore, the fields are of the form. 

u(x,z)~f(z) exp (—ikSx) (30) 

where j(z) is some function of z. The usual assumption in magneto-telluric studies is to set 
S=0 insofar as the fields in the earth are concerned. Cagniard [1953, 1954] justifies the 
assumption by saying that kx is always small compared with unity; consequently, exp (—ikSx) 
may be replaced by unity. The fallacy in his argument is that S may be very large when 
the source of the field is in the near zone. Cagniard does not recognize this important fact 
and, rather surprisingly, he imagines the source field to be a uniform plane wave at a real 
angle of incidence. In this trivial case S is then the sine of the (real) angle of incidence and, 
of course, it could never exceed unity. 

To demonstrate that S may be greater than one, we need only choose a very simple model 
for the source. It is a uniform line current / at some height H above the surface of the earth. 
For simplicity, it will be located at z=—H and runs parallel to the y axis of the coordinate 
system (x,y,z). The primary field (neglecting the influence of the earth) of this line current 
has only y component of the electric field. It is well known and is given, in terms of a modified 
Bessel function, by [Wait, 1959] 

EjJ&f Ko[*B] (31) 

where R=[(z+H) 2J r x 2 ] i is the distance to the current line from the observer at (x,z). If, 
rather hypothetically, &J?>>1, the Bessel function may be approximated by the first term 
of its asymptotic expansion. Then 

E '=-%r(2-ss)'«-"- < 32 > 
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Tims the magnitude of the field is proportional to (kR)~^ and the wave front is nearly plane. 
However, unfortunate as it may be, this limiting form is never achieved in the micropulsation 
region. For example, at co/27r--10~ 2 the wavelength 27r/&^3X10 7 km, which is rather large. 
In nearly all cases of practical interest the observer is in the near zone where kR<^<^\. Thus 

^_j/g/ [log {kRj2) +05772 _ _ >]f (33) 

which bears little similarity to a plane wave. 

To demonstrate the existence of S values greater than unity, it is desirable to write the 
equation in the form of an integral [Wait, 1953b], thus 

E -it^. J* exp [-ikC\z+H\] e _ USxdS {M) 

where C=(l— S 2 )K In the far field, where (&7?)>>1, the important values of S are in the range 
between and 1 and the integral can be evaluated by the method of stationary phase to yield 
eq (32). However, if kR is small the important contributions in the integral are for large 
values of S. 

To focus attention on the physical aspects of the problem, a single harmonic component of 
the spectrum is considered. Also, initially the earth is assumed to be a homogeneous half- 
space of conductivity (r x and permeability ju . The earth curvature is neglected and the assump- 
tion is justified in the appendix. Thus, for a source field which varies as exp (—ikSx) and 
when the electric field has only x component, the surface impedance is easily found to be [Wait, 
1953b, 1958] 

Zi^-EJH&^Q^y (l-fc)-» (35a) 

where 

/ 3=F5 2 /<r 1 M a ) . 
If i3«l 

Z^(^f (35b) 

which is the value appropriate for a homogeneous half-space under the assumption of negligible 
horizontal gradients. Since kS can be identified physically with d/dx it follows that the condi- 
tion j3<C<0 is equivalent to saying that the surface fields vary in the x direction in a distance 
small compared with 8 where 

J=[2/(oi/i«)]* 

is the skin depth of the conducting medium. At a frequency of 0.1 c/s and for a conductivity 
ci^lO -3 mhos/m, 8 is approximately 50 km. Consequently, if the fields vary appreciably in a 
distance of the order of 50 km over the surface of a homogeneous earth, some departure from 
eq (35b) is to be expected at frequencies less than 0.1 c/s. This limitation was pointed out 
previously [Wait, 1954] using a somewhat different argument. 

When the earth becomes horizontally stratified, the limitation that f$ is small remains. 
However, the stringency of this condition varies in a manner depending on the stratification as 
pointed out by Price [1962]. To demonstrate this interesting phenomenon, the surface im- 
pedance for a two-layer earth is considered. The model consists of a homogeneous surface 
layer of conductivity <j x and thickness h } and semi-infinite lower layer of conductivity cr 2 . Then 

on the assumption that ^~ =—ikS and ^-=0, it follows from previous work [Wait, 1953b] that 

Zi=(^)\i-tf)-*e (36) 
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where 

n _ 67+tanh x , 

^~l+£tanhx {6/) 

where 






and 

X=(wimw)^i(1— ip)K (39) 

If Ai tends to infinity, Q approaches unity and eq (36) reduces to (35a). Thus Q can be regarded 
again as a correction factor which accounts for the stratification in the earth's crust. If p is 
replaced by zero, the correction factor becomes 

~ = (o-t/o^+tanh [(WnQ^il ^ / 4Q x 

^ 1 + (<ti/<7 2 )* tanh [(iaiii^h^ 

which is a special case of eq (25) when A 3 =oo . 

To illustrate the influence of finite p on the behavior of Q, a number of calculations were 
carried out on a digital computer. Extensive tabulations of this quantity for a range of values 
of the parameters are available [Jackson et al., 1962]. Some of the results are shown graphically 
in figures 10a to 12b. In the first pair, a 2 is effectively zero corresponding to a conducting 
layer (of thickness h) lying on an insulating substratum. It is apparent from these curves 
that a finite value of p leads to a major change in the shape of the curves. The same behavior 
is evident when the lower layer is finitely conducting as can be evidenced in figures 11a and lib 
where a l fa 2 =25. However, in this case, the effect is not so pronounced. In fact, when the 
lower layer is relatively highly conducting, the influence of finite p is relatively small as can be 
seen from figures 12a and 12b. 

To discuss the significance of the results shown in figures 10a to 12b, it is desirable to 
introduce a scale distance L defined bv 



=^4 



It is a measure of the horizontal distance in which the field changes by an appreciable amount. 
For a periodic disturbance, L is the (spatial) wavelength. Thus 

/3=5 2 /(2Z 2 ) 

where 6 is the skin depth in the upper layer. For a highly insulating substratum, it appears 
that P must be very small or that 5 must be much less than L if the standard magneto-telluric 
interpretation is to be applied. For example, if 

Ai=1.4 km, o-!=10" 3 mhos/m, /=0.1 c/s and o- 2 <<cri, 
it follows that 

6=50 km and (<rijtico)%i=V2Ai/$^0.02. 

Then from figure 10a it can be seen that if j8<10 -3 , \Q\ is within 10 percent of its value for 
P=0. This condition is equivalent to 

L>8/(P^[2) 
or 

L>l,100km. 

The condition is even more stringent for the phase. 

Uniformity of the exciting fields over distances of the order of 1,000 km would not be very 
common, particularly in higher latitudes where the currents could be quite localized. Also, at 
low latitudes the presence of equatorial electrojet [Vestine, 1960] at heights of the order of 100 
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km would be expected to produce considerable nonuniformity. Actually, the two-layer model 
with a highly insulating substratum is a most unfavorable circumstance. Furthermore, it is 
somewhat hypothetical since the poorly conducting crystalline rock will also be of finite depth. 
Consequently, a three-layer model should again be adopted for a more realistic appraisal. 
In this case, eq (36) may be generalized [Wait, 1953b] to 

A 

n _ ffQ+tanh X 

*% A 



1+GQ tanhx 

where 

a A 

A_ (?+tanh x 

1+G tanhx 

gJ 't 1 "^ and X=(i<T2w)W-ifai/*2)K 

L "3/ a 'l~"^Pj 

Now, if the intervening stratum is poorly conducting such that 

(o- 2 Mow)^2«l and /3VowoA 2 «l 
it readily follows that 

Q _ GG+ttmhx 

~l+&? tanhx 
where 



w- ij- 



L<r*/*i—iP- 



The latter equation for Q does not depend on <r 2 or A 2 and thus the three-layer structure is 
equivalent to a two-layer model whose constants are <r u c 3 , and hi. If, m addition, /?<<0 this 
formula for Q reduces to eq (29). 

In summary, the restrictions on the use of the curves of Q given in figures 2a to 9b can be 
applied to a three-layer structure provided the single condition 

i^ 2 i«i 

is met. Thus, for highly insulating substrata, the condition can become quite stringent since 
\Q 2 \ is then large compared with unity. On the other hand, for highly conducting substrata, 
\Q 2 \ may be small and the condition is not at all stringent. 

5. Appendix. Surface Impedance of a Spherically Stratified Conductor 

In many boundary-value problems involving waves in stratified media the solutions may 
be quickly obtained if the analogies with transmission line theory are exploited. In general, 
the analogous transmission line is nonuniform in the sense that the characteristic impedances 
may be different in the two directions on the line. Although Schelkunoff [1943] has discussed 
nonuniform transmission line theory, it seems worth while to include a short exposition of its 
essential features. 

The starting point is the equations which connect the transverse voltage V between two 
parallel wires and the longitudinal current / in the lower wire. In terms of distance x down 
the line, these are 

^=-ZI (Al) and^=- YI (A2) 
dx dx 

where Z is the distributed series impedance per unit length and Y is the distributed series 
admittance per unit length. Allowing Z and Y to be both functions of x one readily finds 
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that 

d 2 V Z' dV TZVTr n , AoN 

1^~Y Tx~ YZV=0 - ( A3) 

Second-order differential equations of this type possess two linearly independent solutions. 
The general solution is a linear combination of these solutions. For example 

V(x)=AV+{x) + BV-(x) (A5) 

where A and B are independent of x and V + and V~ are the fundamental wave functions. 
Similarly, 

I(x) =AI+ 0) + BI- (x) (A6) 

where I + and /" are the corresponding fundamental current wave functions. It is evident 

that the wave functions individually satisfy eqs (Al) and (A2). 

A characteristic wave impedance may now be associated with each pair of wave functions; 

thus 

V + (r) 1 (II + ZV+ 

K+(x) = = — — = — (A7) 

w I+(x) YI+ dx dV+/dx l ; 

and 

K ~( x) =~I=ld = YF lx' == dV ir /dx t (A8) 

It is a convenient physical artifice to consider the waves associated with K + as propagating 
in the positive x direction and those associated with Kr as waves propagating in the negative 
x direction. In the case of a uniform line where Z and Y are constant, the meaning of these 
wave functions is quite clear. In the presence of any local nonuniformity, a reflected wave 
would be generated and the individual wave functions no longer are purely propagating. 
Nevertheless, in many cases, the elementary wave functions bear considerable resemblance 
to "traveling" waves. For example, this happens when Z and Fare slowly varying functions 
of x. 

For present purposes the label progressive is used to describe a fundamental wave on a non- 
uniform line. In the limiting case of a uniform line the progressive waves become traveling 
waves. 

It is important to note that there is some arbitrariness in the selection of the wave functions 
on nonuniform lines. The choice is usually made on the basis of convenience. 

A useful quantity is the ratio of the wave functions at two points x x and x 2 on the line. 
Thus, by definition, 

Xv(X 1 ,X 2 )=y+, y Xv(Xl,X 2 )=y_, v (A9) 

Xi{x u X 2 )J^y *?(*!, Z 2 )=j^g- (A10) 

We are now in the position to study the reflection in nonuniform lines. For example, a 
semi-infinite nonuniform line is terminated in an impedance Z. If Z=K + , the voltage associ- 
ated with this incident wave is exactly equal to the voltage across the terminal impedance and 
thus the entire incident current wave flows through this impedance. Hence, no reflection 
occurs. However, if Z^K + , the incident current wave is not completely absorbed and reflec- 
tion occurs. If the incident "progressive" waves are characterized by the wave functions 
V + {x) and I + (x), the problem is to calculate the reflection coefficient at the terminal impedance. 
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These are defined by 

J- v 

2/=j+ and Q.v=y+ (All) 

for current and voltage, respectively. Now, 

V t =V++V~ and I t =I++I~ (A12) 

where V t and I t are the voltage current at the terminals of Z t . Furthermore, 

V+=K+I+, V~=K-I~, V t =ZJ t . (A13) 

The latter two sets of equations are readily solved to give 

M+—Y 1 1 

qv= M~+Y where M±== K± and Yt= Z t ^ Al4 ^ 

and 

K + -Z t , A1 _. 

^ir+z; (A15) 

A basic problem in nonuniform lines is to form wave functions V(x) and I(x) at x=Xi in 
terms of the specified impedance at x=x 2 . For convenience, x 2 ^>x. Using the preceding 
conventions one may easily write 

V(x)=VHx)+V+(x 2 )q v (x 2 )¥^ ( A1 6 ) 

and 

I(x)=I+(x)+mx 2 ) qi (x 2 )j^y (A17) 

Here V + and I + are regarded as the "incident" progressive wave and V~ and I" are the "re- 
flected" progressive wave. 

The impedance at the point x=x x is then given by 

?(„ \— y( X V — TS+f™ \ l + gF(^2)X^(Xi,X 2 )Xy(x 2 ,Xi) /i 1Q \ 

A{Xi) -JM- K Wl+toixJxtiXirixtiw)' (A18) 

This result is most useful in reflection-type problems. 

To complete this very brief survey of nonuniform transmission line theory, the correspond- 
ing formulas for transmission coefficients are also given. 

The transmission coefficients at impedance discontinuity are defined by 

Pv=yi and Vi=j± (A19) 

being analogous to the previous definitions of q v and q T . Then from eqs (A 12) and (A13), 



and 



^M-+F7 (A20) 



»-z=+zr- (A21) 



As expected, these coefficients become unity when Z t =K + (or Y t =M + ). 

The transmission coefficient across a section (x h x 2 ) of another line inserted between the 
original line and impedance Z t is readily obtained. The result, given by Schelkunoff [1943], is 
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T=p(l + qX+qW + . . . )x+(x ly x 2 ) 

V-X+^Xt) (A22) 



1 — qx 
where 

p=p + (x 1 )p + (x 2 ), q=q-(x{)q + (x 2 ) (A23) 

and X=X + (xi,£ 2 )x~( ;r 2, ;r i). The result holds for both current and voltage waves which ex- 
plains the absence of subscripts on / and V. The physical meaning of the results is very 
clear: p is the product of the two transmission coefficients at the discontinuities x x and x 2 , g is 
the product of the reflection coefficient for a progressive wave incident from the right on the 
junction x u and the reflection coefficient for a progressive wave incident from the left. (Here 
x increases towards the right.) In the absence of multiple reflections T would be equated to 
px + (xuX 2 ). Thus, the factor (1+#X+<Z 2 X 2 + . . .) represents the influence of multiple 
reflections. 

We now proceed to apply nonuniform transmission line theory to a spherically stratified 
medium. In particular, let us consider the earth (of outer radius a x ) as consisting of a homo- 
geneous core of radius a 2 of electrical constants a u e h and ii\ surrounded by a homogeneous man- 
tle of electrical constants a 2 , *2, and n 2 . It is assumed that the sources of the field are completely 
exterior to the earth. Thus, within the concentric homogeneous regions the fields may be 
derived from two scalar potential functions U and V. The first set which are TM (transverse 
magnetic) are derivable from U and the second set are TE (transverse electric), derivable 
from V. 

For TM waves 

while for the TE waves 



Both U and V satisfy 



E ^h(^ u )' H <=° (A24) 

* U ' 1 *Uem+'™=^U (A26) 



dr sin 6 dd 



as can be readily ascertained from Maxwell's equations. 

Since any field may be expressed as a superposition of these two sets, it is sufficient to 
discuss them separately. 

Using the standard separation-of -variables technique, it is assumed that 

U=u(e#)u(r). (A27) 

On substituting this into eq (A26) one finds that u and u must satisfy 



and 



sin 6 ~ (sin ^+^+,(,+ 1) S in 2 du=0 (A28) 

ff-[T 2 +^]^0 (A29) 

where v is a constant. 

The u's are Legendre functions and the general form may be expressed as 

u(0,*) = S ([Fm(v)P? (cos o)+G n (v)P?(-cos 6)] dvX[F m cos m4>+g m sin w4] (A30) 

m J 

where the integration contour is suitably chosen in the complex v plane and m is an integer. 
The constants F m (v), G m (v),f m , and g m depend on the nature of the source field. For example, 

535 



if the source is a vertical electric dipole at 6 = 0, F m (v) = Q } and g m =0. Also, for m^O, / m =0. 
Tims 

u= Cg (p)P„(-cos $)dv. (A31) 

It is known that integrals of this type can be deformed so as to enclose the poles of 6o(v) and 
thus u can be represented as a series of residues. Therefore, 

us*2*i X! [Residues of G (v 8 )] P„ s (-cos d) (A32) 

Va 

where v s are the solutions of [G (v)]~ 1 = 0. Provided |^ S |>>1 and 6 is not near or -k 

P v (—cos 6) ^ c . on ^ \ e-ika.es^ei^e-tka^^-eisl 1 (A33) 

sK J (sin 6)i L J 

where v s (v s -\-l)^v*^ka,iS s . Here the angular function u has the physical character of two 
waves traveling in opposite directions around the cylinder. In analogy to waves on a flat 
surface, S s can be interpreted as the sine of a complex angle of incidence where k is the wave 
number in free space. 

It is clearly apparent from this simple example that v is related to the azimuthal variation 
of the source field. In the general case, v depends both on the longitudinal and latitudinal 
variations. 

Solutions of the radial equation are conveniently expressed in the form 

u=AI v (yr)+BK v (yr) (A34) 

A A 

where I v and K v could be any two independent solutions of eq (A29). Herey=[inu(<r-{-ieoo)Y where 

A 

Re 7> 0. For convenience I v (yr) is chosen to be a solution which is finite or zero at r=0 and 

A 

K v (yr) is chosen to vanish as r tends to infinity. In terms of modified (cylindrical) Bessel 
functions 



and 



^(«)=(^)*/, +1/1 (a) (A35) 

K v (z)=(^J K v+1/2 (z) (A36) 



where / and K have their conventional meaning. 

In the present problem, the waves associated with I v (yr) are regarded as being "incident," 

A 

whereas the waves associated with K v (yr) are regarded as reflected." 
The surface impedance at r=a x is defined by 

^=f-S I (A37) 



L J^ V (f>Jr = a 1 



which, in general, depends on v. The medium between the limits r=a x and a 2 is now regarded 
as a nonuniform transmission line of length l=a x —a 2 . The transverse voltage V T and the 
current / on the line are then analogous to the electric field — E v e and the magnetic field H^, 
respectively, where the superscript v is to indicate the possible dependence on v. 
The characteristic impedance of the line looking inward is then 

^ + (7ir) = r7iP^ (A38) 

A(7ir) 
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and the impedance looking outward is 



K~(y 1 r) = -vi 



A 

K v (7ir) 



(A39) 



In the above 
and 



7i=['lAtiw(o'i+'i€ico)p 

The line is now considered to be terminated by an impedance Z t where 

/i(Tir)" 



Z t = 



V2; 



L Uy?). 

From the analogy with transmission line theory, one readily finds that 

Z =K + ( a ) [l + qpUa2,ai)x e (a u a 2 )^] 



where 



-q fl X h (a 2 ,a 1 )x h (a l ,a 2 )_ 



(l/Z t )-l/Jg+( 7 ifli) 
qe (l/Z t ) + llK-(y iai y 

= Z t -K+{ yi a x ) 
qh Zt+K-iy^)' 



x e {a 2 ,a l )-- 



dJ v {y Y a 2 ) 
ajyiyidx) 



x e (a h a 2 ) = 



a 2 K' v (yi(ii) 



aiK' v (x\a*) 



x,(a 2) a 1 )=^^- ) , 

^2/,(7l«l) 



x h {a\,a 2 )-- 



_ a 2 K v (y l ai) 

A 

diKyiyfa) 



The preceding results can be greatly simplified under the assumption that y^ 

A 

are large compared with unity. For example, noting that I v (z) satisfies the equation 



d 2 t(z) 
dz 2 



(i+^ti)]/, ( ,), 



(A40) 

(A41) 
(A42) 
(A43) 

(A44) 
(A45) 
(A46) 

(A47) 
and y 2 a 2 

(A48) 



it readily follows that K + (yr) satisfies 

*,+„ f =[,+*+!)],;, (A49) 

where K=K + (y t r) and z=y x r. For a first approximation, the derivative term may be neg- 
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lee ted. Thus 



* + <->«4 i+ wT 



(A50) 



This approximate result, obtained in a very simple fashion, corresponds to the use of the Debye 
or the second order result. The theory of Bessel functions indicates that its validity, in this 
form, depends on the conditions Re 7ir^>^>l and |7i?*| 2 somewhat greater than v(v-{-l) [Sommer- 
feld, 1949]. Similarly, 

^( 7l r)^^[l+^^J (A51) 

and 



K ^-^b+Wl 



(A52) 



To within the same approximation, the reflection coefficients at r=a t may then be written 

„_ L_\72«2 /J L \7l&2 /_ 



and 



(XM 



<i„ 



\T2a2 /J L \7ia 2 / 



(A53) 



*72 


_ \72a2 / _ 


2 
—771 


\Tl«2 / _ 


V2 


- \y 2(12 J _ 


2 





(A54) 



where v(v+l)^v 2 ^— y 2 S 2 =k 2 S 2 . The corresponding (approximate) form for the product of 
the transmission coefficients is easily shown to be 



X e (a2,a 1 )x e (a 1 ,a 2 )^x h (a2,ai)x h (a 1) a2) 



(A55) 



where l=a x — a 2 . 

Using eqs (A50), (A53), (A54), and (A55), eq (A41) for Z v can be expressed in a fairly 
convenient form despite its complicated appearance. A great simplification can be made 
when the thickness I of the shell is small compared with a x . Thus, the ratios aja 2 and a,i/r 
in the preceding expressions may be replaced by unity. The surface impedance may then 
be written 

|[|^ST+tanh [(yl-ytS 2 )H] 



**-1i\} 7? ^J .y\ [yl-ylS 2 y 
7^[_Y?-7^ 2 J 



tanh [{y\-ylS 2 )H} 



(A56) 



where it has been assumed that Mi=/*2=Ato- This result is identical to what one would obtain 
on a planar two-layer earth for a plane wave incident at an angle of incidence arc sin S [Wait, 
1958]. Consequently, in this limiting case the influence of earth curvature vanishes. 

The specific results developed above for the surface impedance are restricted to TM 
waves. The corresponding results for TE waves are obtained by exactly the same method. 
In fact the results are completely analogous if admittances are used in place of impedances. 
The surface admittance defined by 

V — — 
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has precisely the same form as eq (A41) for Z v if K + {y x a) is now replaced by the admittance 
M + (yid) and q e and q h have the same form as eqs (A53) and (A54) except that ^ and 772 are 
replaced by their reciprocals everywhere. The bracketed terms are unchanged. 

Under the assumptions that Ke 7^1 and Re 7!a 2 >>l, Mi=M2 = Mo, and that a x — a 2 <<ai, 
one easily finds that 

" 7 !+[^g] tanh[ W -^)«] 

This is equivalent to the result given by eq (36) in the main body of the text. 
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